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I. INTRODUCTION 



Electromagnetic polarizabilities are fundamental properties of hadrons which are man- 
ifest in various processes involving two photons Recent measurements of the 
electric, a, and magnetic, /3, polarizabilities of the nucleon narrowed considerably experi- 
mental uncertainties in these observables, and were accompanied by a number of theoretical 
studies. Attempts to describe a and f3 were made in various approaches, ranging from chiral 
perturbation theory and dispersion relations M to chiral soliton models [p|-|i4|j. Earlier 
calculations of a and (3 are reviewed in Ref. |fT5[| . 

In this paper we calculate the average electric polarizability of the proton and neutron, 
a = \{a p + a n ), using the solitonic approach [0,0 to the Nambu-Jona-Lasinio model |L8 



This model has led to quite successful phenomenology p9| , p0| . In particular, baryonic mass 
relations, magnetic moments, and various form factors are reproduced well within expecta- 
tions. It is therefore challenging to see if we can also describe two-current observables, such 
as the polarizability, in the NJL model. 

There are two important features in our study which have not been considered in earlier 
works in soliton models: 1) the role of the Dirac sea effects (this can only be done in a model 
which has the Dirac sea, such as the NJL model), and 2) inclusion of rotational l/iV c -effects. 



Such effects were recently analyzed for the case of the axial coupling constant g& in Ref. f21 
and were more closely examined in Refs. p2lp3|. 



In the NJL model the only dynamical degrees of freedom are quarks, which occupy 
valence orbitals, as well as the Dirac sea. We find that the effects of the Dirac sea in our 
calculation of a are very well reproduced by the first two terms in the gradient expansion, 
where the first term is the so called "seagull" contribution to a, discussed in many previous 
papers |9,1C,12]. Our analysis demonstrates explicitly that the inclusion of the sea-gull in 
the cr-model or in the Skyrme model does not violate gauge invariance, as has recently been 
claimed in Ref. [23J]. This is readily seen from the full NJL expression, which is purely 
dispersive in character (i.e. involves no seagull terms) and manifestly gauge invariant. The 
seagull term emerges from the full expression when the gradient expansion is performed (see 



Sec. pj). 

The inclusion of rotational 1/N C effects has very sound phenomenological consequences. 
As discussed extensively in Ref. [[T2|], chiral soliton models, when treated at the leading- N c 
level, have problems in predicting electromagnetic polarizabilities correctly. The dominant 
contribution to the electric polarizability is obtained from pion tail effects, and is propor- 
tional to g\. If one insists that the model fits the value of g& (which one should!), then 
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the value of a obtained at the leading N c level is roughly a factor of 2 — 3 too large, as 



demonstrated in the cr-model calculation of Ref. [ I2| . In the NJL model (as well as in the 
Skyrmion) the value of obtained with various fits has always been notoriously too small 
at the leading-iVc level ||25|| , i.e. of the order 0.7 — 0.8. As first noticed in Ref. [pip , the 
inclusion of rotational 1/N C corrections in the NJL model is capable of raising qa to a com- 
fortable value of ~ 1.2 [ p2"p3| . This is a big correction: it raises g\ by a factor of 2 — 3. 
We calculate analogous 1/N C corrections for the electric polarizability, and find a sizeable 
contribution. As a result, and after including approximately additional corrections due to 



the N-A mass splitting [2q|, we are able to obtain a number which is closer to experiment 



than in other studies in soliton models fl2| , |14|1 , but still too large. For the typical choice 
|T§| , |2T)|| of the constituent quark mass M = 420 MeV we obtain a ~ 19 x 10~ 4 fm 3 compared 
to the experimental value a exp = 9.6 ± 1.8 ± 2.2 x 10~ 4 fm 3 f|-|5]]. 

In this paper we do not analyze the splitting of the proton and neutron polarizabilities, 
since it involves a complicated problem of treating the translational zero mode Also, 



the magnetic polarizability is not analyzed. It was shown in Ref. |12] that the large-iV c 
approach is not a good starting point to describe this observable. Hence, we concentrate 
solely on the average proton and neutron electric polarizability, a. 

The outline of the paper is as follows: In Sec. [H] we develop the necessary formalism 
of linear response theory for the NJL model, and derive a basic expression for the electric 
polarizability. In Sec. |TJ we calculate the leading- A c contribution to a, which comes from 
both valence and sea quark effects. In Sec. [TV] we calculate the rotational 1/N C corrections 
to a. The valence contribution is calculated exactly, while the sea part is estimated using 
gradient expansion. Additional corrections, due to A-N mass splitting, are discussed in 
Sec. |V[ Sec. [VT] contains our results and conclusions. 

II. FORMALISM 

Polarizabilities are defined as coefficients in the low-energy expansion of the Compton 
amplitude It has been shown |l],[|,|24| that the polarizability measured in Compton 

scattering, a, can be written as a = a + Aa, where 



„ | (N | D z | X) I' eV) E eV + 1) 

a = 2 £ N En — Ex ' = ~3M W 



Here \N) is the nucleon state, and \X) is any intermediate excited state connected by the 
electric dipole operator D z . The first term in the expression for Aa is the so called recoil 
correction, which involves the charge, e, the mass, M, and the mean squared charge radius 
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of the particle, and the second term is the Schwinger scattering term for a particle with 
anomalous magnetic moment k. Our goal is to calculate the dispersive contribution a in 
Eq. (TJ). Note that according to Eq. ([!]), this is equivalent to calculating the coefficient of 
the second-order energy shift of the nucleon in the constant electric field. 



First, we very briefly review the NJL model |16| , ^9| . Our expressions are written in the 
Euclidean space-time. In the leading- iV c treatment the effective action results from taking 
into account one-quark loops. It contains the quark part, S^ s , the chirally invariant mesonic 
part, S^g, and the chiral symmetry breaking part, S^, which provides the pion its mass 

Seff = S^g + S^fi + S^g (2) 



S*s = -Tr ln A (^ - Mf/ 75 ) - vac, (3) 
J d 4 x (a 2 + Tr 2 ) - vac, (4) 



a 2 

qM _ P 
D eG — TT 



b f^xa-vac. (5) 

9 J 

The trace is over the quark spinors. The subscript A reminds of the NJL cut-off (in this 
paper we use the proper time regulator p7|j28| ), m is the current quark mass, and M is 



the constituent quark mass. The matrix U Tj represents the meson hedgehog profile in the 
non-linear representation 

W 6 {r) = exp(i7 5 r ■ fO(r)) = —(a h (r) + i-f 5 r ■ 7Th(r)), (6) 

where O(r) is the chiral angle, o~h(r) = crft(r) and 7T/j(r) = rir^r) are the sigma and 
pion hedgehog profiles with 07,, (r) = f n cos O(r) and 7r^(r) = f n sin O(r), r is the radial unit 
vector, and r are the Pauli isospin matrices. The parameters of the model are fixed by 
reproducing the experimental values of the pion mass, m n , and the pion decay constant, f^, 
leaving the quark-meson coupling constant g, or, equivalently, the constituent quark mass 



M = gf n as the only free parameter [19 



Solitonic solutions to the NJL model are found by selfconsistently solving the coupled 
Euler-Lagrange equations for the profile function O and the Dirac spinors q n , which satisfy 
the Dirac equation 

Hq n = (-iV • a + (3MU 15 ) q n = e n q n (7) 

In the presence of the the electromagnetic field the minimal substitution 

_|_ ieQA 11 , with the quark charge Q = ( + |t 3 ), leads to the following gauged 
effective action 
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= -Tr ]n A (i# - eQ4 - Mf/ 75 ) - vac. (8) 

The hedgehog solution brakes the spin, J, and isospin, /, symmetries of the lagrangian, 
preserving the grand spin G = J + I symmetry. To restore good quantum numbers we use 
the semiclassical "cranking" projection scheme p9| , |30|| . Introducing the collective rotation 
matrix R(t) = B exp(i^^-t) in the usual way and transforming the quark fields and meson 



profiles to the isorotating frame [^,0 we obtain 

= -Tr ]n A (i0 - eB ] QBfi - ~ MU 1B ) - vac. (9) 

It is convenient to introduce the notation B^t^B — c-t, where c a is an element of a Wigner 
D-matrix, c a = _D* 3 [|3T|. Then the quark electric charge in the isorotating frame can be 
written as 

The action (§) describes the system which is perturbed by small perturbations: electro- 
magnetic, and cranking. From now on our method will follow very closely Ref. [31]. The 
idea is the following: we have found the soliton, and now we want to calculate its linear 
response to external perturbations. The perturbation operators consist of products of two 
parts: intrinsic, which act on the internal structure of the state, and collective, which act 
on collective coordinates. For the cranking perturbation, iflQ^-, the collective operator is 
f2, and z/3r/2 is the intrinsic operator. Similarly, the electromagnetic perturbation has a 
collective part (vector c in Eq. (JlQ|)) and an intrinsic part. Also, the nucleon state is a prod- 
uct of a collective and an intrinsic part [|29| |3Tj] . Formal perturbative expansion of Eq. @ 
to n-th order in small parameters involves expressions of the general form 



«ix,..,<n = (coll\V™ ll V% ll ...V™ ll \coll) 

J d 4 Xl ... d 4 x n (H\V^ tr (x 1 )G(x 1 ,x 2 )V^ tr (x 2 ) 



G(x 2 ,x 3 ) ... V^ 1 (a;„_ 1 )G'(a; n _ 1 ,x n )Vr r (x n )G'(x n ,a; 1 )|iJ), (11) 

where \coll) is the baryon collective wave function, \H) is the intrinsic hedgehog state, Vi k is 
the perturbation of the k th type, and G is the quark Feynman propagator in the background 
meson field, which in the spectral representation has the form 

G(x,y)= (e(x 4 -y 4 ) £ e-^-w) - 9(y 4 - x 4 ) £ e—^) ) qn ( x ) q l( y ), (12) 

\ nSunocc. ragocc. / 

where occ. denotes all occupied states, i.e. the valence as well as the Dirac sea states and 
unocc. denotes the unoccupied positive energy states. 



Examples of application of Eq. ( |TTD are given in Ref. [31 
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III. LEADING- N c CONTRIBUTIONS TO a 



It was shown in Ref. |T2[ that the leading contribution to the electric polarizability of 
the nucleon scales as N c . Our first task is to obtain this leading contribution for a in the 
NJL model. We denote it by a^°\ Since the cranking frequency fi ~ 1/A^, the leading- N c 
contribution is obtained by expanding the effective action (§) to second order in the electric 
field, and to zero order in Q. From Eq. (|T0|) we notice, that for the leading N c contribution 
we have to keep only the isovector part of the charge, since the isoscalar part is one power 
of N c down. 

We now take A 4 = Ez, A = 0, and expand Eq. (§) to second order in E and zero order 
in O. We obtain the formula of the form (|Tl[) 



a (0) = __ e 2 {N | Tr (^JL X3 Q ^ y) £_I G{ y jX) ^ | N) 

1 9 , , „ > . , ^ (m\r a z\n) (nW b z\m) 

-N c e 2 (N\c a c b \N) — U ^ — . ( 13 ) 
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nGunocc. 



where \N) is the collective nucleon state, T is a large time interval, over which the time 
integrations are carried, and \n) and e n denote the eigenstates and eigenvalues of the Dirac 
hamiltonian ([!]). We can now split this formula into the valence and sea parts. The sea 
contribution carries the cut-off function, and the vacuum contribution has to be subtracted. 
Finally, using (N \ c a c b \ N) = ^5 ab , we get the following formulas for the leading N c valence 
and sea contributions to the electric polarizability of the nucleon 

N c 2 ^ (val\r a z\n) (n\r a z\val) 



(o) 

"vai = — e 



E - ' ' ' 7 ; (i4) 



(■> n ^ val e n e va i 



«iea = ^e 2 E R(e m ,e n , A) (m\r a z\n) (n\r a z\m) - vac, (15) 



where R(e m , e n , A) is the proper time regularization function 
#(e m , e n , A) = — — —= [ 



2 2 2 2 

1 f°° ds e~ S6m — e~ S6n e m e~ S€m + e n se~ S6n 

+ s- 



40? J a- 2 s 3 / 2 



(16) 



with A denoting the (proper time) cut-off parameter. The value of A is fit for any M to 
reproduce the experimental value of the pion coupling constant f n using the amplitude of 



the weak pion decay [jnj. The subtracted vacuum contribution in Eq. (15) has the same 
structure as the first term, but the free quark states and eigenvalues are used. It may be 
immediately verified that both the valence and sea contributions to c^ -* scale as N c [12 |. 
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In the NJL model, predictions depending on the cut-off are somewhat ambiguous. As 
we shall see, however, in our case is dominated by the lowest term in the gradient 
expansion, which is independent of the cut-off. 

The numerical results for the valence, a^j, and Dirac sea, a^l, parts of are plotted 
in Fig. 1 as a function of the constituent quark mass M. Both contributions decrease with 
increasing M, the valence curve has much larger slope for smaller M due to the fast increasing 
energy gap between the valence orbital and the positive continuum. For the optimum value 
T|,|I1 of M = 420 MeV, we obtain a (0) = 16 x 1(T 4 fm 3 . The sea contributes to ct (0) 65%, 



and dominates over the valence part, which gives 35%. Since the Dirac sea effects describe 
the physics of the pion cloud, this result is in qualitative agreement with the results of other 
models [0,0,0. 

We have also performed the gradient expansion of aj® up to first two terms following the 



scheme of Refs. [p3H 35|] . This provides a link to the a-model |3q , |i2| as well as an independent 
check for our numerical procedure in the case of large size solitons. The leading-order term 
in the gradient expansion of gives us exactly the seagull term of the a-model 



"seagull = — ^— J drr A {it h {r)) . (17) 

Note that this seagull term arises in the NJL model from purely dispersive Dirac sea effects 
(cf. Eq. (0)). Recently, the inclusion of the seagull term in soliton models such as the 
cx-model or the Skyrmion has been questioned The NJL model provides a simple 



interpretation of how Eq. (|T7|) arises in the a-model, and clears the doubts of Ref. [24 
Numerically the seagull term gives more than 80% of the total sea contribution. 

Looking at absolute numbers it may seem striking that our value of the seagull contribu- 
tion, ~ 9 x 10 -4 fm 3 , and the <r-model [12] value, ~ 30 x 10 -4 fm 3 differ so drastically. The 
explanation lies in the large difference in the values of the axial coupling constant at the 
leading N c level. Our model predicts g^ ~ 0.75, while in the cr-model g^ = 1.42 [^Oj. The 
asymptotic pion tail, which gives the dominant contribution in fll~7|), is proportional to g^ 



and since the pion field enters quadratically in (|T7j) the large effect follows. We will come 
back to the question of in the next section. 

It is interesting to investigate in some greater detail how well the gradient expansion 
works for the electric polarizability. In Fig. 2 we compare the results of the gradient ex- 
pansion with the full calculation of for large-size solitons parameterized by the profile 
function involving a size parameter R 

0(r) = -2atanf^-(l + m 7r r)e- m ' rr y (18) 
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With increasing size of the soliton the result of the gradient expansion for a™ up to first 
two terms converges to the full result very well. At R — 2 fm the difference is only 2%. 
Even the leading order of the gradient expansion, a seagu n, shows good agreement with the 
full result, reproducing it up to 90% at R = 2 fm. 



IV. ROTATIONAL CORRECTIONS TO a 



Now we turn to the problem of the next-to- leading order corrections in the 1/N C expan- 
sion of the electric polarizability a. As already mentioned above, the value of in chiral 
soliton models is crucial for the results for a. The prediction for in the NJL model (as 
well as in the Skyrme model) at the leading- N c level, denoted by g^, turns out to be much 
too small, ~ 0.7 — 0.8 f25f, and no parameter fits or extensions of the model were able to 



raise this value. Recently, Ref. |21| proposed how to include rotational 1/N C corrections 
to gA in the NJL model. A correct and detailed study of this effect has been presented in 
Refs. p2]J2B|J57| . The nature of these corrections is very easy to understand, for instance 
from the point of view of the formalism of Ref. |3l| . At the leading- N c level the expression 
for g^ involves a sum of matrix elements of the axial operator in the quark single-particle 



states |25|. However, in the presence of cranking the hedgehog state is perturbed, and we 
have a situation where one of the perturbations is the axial vector current, and the other 
one is the cranking. The effect discussed in Refs. pl,23,22| is associated with the noncom- 
mutativity of the collective parts of the corresponding operators (the cranking and the axial 
vector current operators in this case). An analogous effect arises for the case of a, as we 
discuss below. 

Since the cranking operator scales as 1/N C (it involves the cranking frequency ft ~ 
the correction is suppressed by one power of N c . This correction to the axial vector coupling 
constant, denoted by g^', has the valence part and the sea part. The numerical value for g A ^ 
is about 0.5 pl|-f23|, and raises the total value of gA to a comfortable number. Of course, the 
described 1/N C corrections are not the only possible ones, but they seem to be particularly 
important. 

We have explained before that the small value of reduces the value of the seagull 
contribution to the polarizability, which is a desired effect. Since roatational 1/N C corrections 
raise significantly the value of gA, we have also to examine how analogous effects influence 
the electric polarizability. The 1/N C rotational corrections to a are obtained by expanding 
the effective action @ to second order in E and first order in O. Thus, we have formally 
the case of the third-order perturbation theory. We get expression of the generic form (|TT|), 
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with all possible time-orderings of collective operators. 

The relevant collective matrix elements for the present case have the form 

(N | c a Q b c c | N) = —e abc , (19) 

where J is the moment of inertia. Finally, after straightforward algebra, we obtain for the 
valence contribution 

iN e 2 e abc 



ft 



(1) ._ 

val U j 



(20) 



X 



yv (val\r a z\n)(n\r b \m)(m\r c z\yal) y-, (val|r a z|n)(n|r 6 |m)(m|r c ^|val) 

- e va i) (e m ^val) n<0 (^n m ^val ) 

m>0,m^val m<0 

^ ( v a-l|T a z\n) (n\r b z\m) (m|r c |val) 2 ( va l| ra £|w)(H r6 ^l m )( m l rC l va l) 

n>o,n^val ( e n e val) ( e n e m) n <0 ( e n e val)( e n e m) 

m<0 m>0 

and for the sea contribution 

i N p 2 afec 

12J 1 J 



(fc|r a z|n)(n|T 6 |m)(m|T c z|/c) ^ y-, (/c|r a z|n)(n|r b |m)(m|r c z|A;) 
^ (ffc — £m)(tk — € n ) k>0 (ek — e m )(efc — e n ) 



k<0 

m,n>0 m, n<0 



vac. 



In the above expressions notation n < means the level with e n < 0, etc. Recalling that 
J7" ~ N c , it can be verified that is of order N®, while the leading contribution to ft^ ^ 
was of order N c . 

The numerical value for ft^j is obtained directly from Eq. (|20|) . It involves a double sum- 
mation over quark levels, and computer time involved is similar to the calculation of, say, the 
moment of inertia. Our results are shown in Fig. 1. The value of a^l is about 3 x 10~ 4 fm 3 
for M = 420 MeV, so the effect is about 60% of the leading- N c valence contribution a^. 
The decrease of a£] with increasing M is similar to the behavior of the leading- N c valence 
contribution a^J. 

Now we turn to the sea contribution oj^. This contribution in principle may involve a 
cut-off, similar to other observables calculated in the NJL model. However, since Eq. Q2"T| ) 
involves a computer-time-consuming triple sum over quark states, we approximate a^H by 
the leading term in the gradient expansion. As a bonus, we get the nice feature that the 
result is independent of the cut-off, exactly as in the case of the leading- N c contribution 
CKseaguii- In the present case the gradient expansion technique involves some tricks p5| , so 
we sketch the method in some greater detail. The first step is to rewrite Eq. fl21~D in a form 
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involving unconstrained sums, i.e. sums over all quark levels. This can be achieved using 
the identity true for any real numbers a, b and c 

du>iduj 2 duj 3 111 1 1 

(27r) 3 w\ — ia lji — ib u)\ — ic iO\ — uj 2 — if] uj\ — uj 3 — ir\ 

= -1^-^ — *— e(-o) e(b) 6(c). (22) 

a — b a — c 

With help of Eq. (p2|) we can replace the constrained sums in Eq. ( f21~D by sums over all 
indices. For instance 

yv (k\T a z\n)(n\T b \m)(m\T c z\k) 



k<0 
m,n>0 



I 

k,m,n 



du)\(ko 2 dia 3 



^2 (k\T a z\n)(n\r b \m)(m\T c z\k) J 



(27r) 3 <jj x - ie k Lo x - ie m lo x - ie n 
P + inS^x - u 2 )) (p + indium - u 3 j) , (23) 

UJl — UJ 2 J \ UJ\ — UJ 3 / 

where P denotes the principal value integral. Introducing the trace over the quark spinors 
Tr we can now rewrite Eq. ([H]) as 

(1) r du,du 2 1 

12J J (2tt) 2 uj x -uj 2 1 J 



xTr 



r z r t z V r z r r z- 



— vac. 



c<Ji — iH uj 2 — iif c^i — iH uj x — iiJ uj x — iH u 2 — iH _ 

Standard gradient expansion techniques may now be used to expand the functional trace 
in Eq. (|4j) J£|. The lowest order term (denoted by i OW est), contains no gradients, and we 
obtain 

(1) 32N C M 2 r d 3 k r doj Y du 2 1 (2ujx+oj 2 ) 

«sea,lowest ~ ^jJT J J^y J ( 2vr )2 ^ _ ^ + f,2 + ^2)2(^2 + fc 2 + M 2) 

x | d 3 r z 2 (7r,(r)) 2 

- ^ M ^ 

" TexTif ttscasu11 ' ( j 

where a seaguU is given in Eq. ([17]). 

For M = 420 MeV the moment of inertia has the value J = 1.1 fm, and we obtain for the 
proportionality factor in the last equality in Eq. (p5|) the value 0.54. The numerical results 
for lowest are plotted in Fig. 1. The contribution increases with increasing constituent 
quark mass M. This is due to the decrease of the moment of inertia with M. Guided by 
our experience from gradient-expanding the leading- N c contribution (Sec. p| ), we hope that 
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also in the present case the lowest term in the gradient expansion almost saturates the value 

O f "sea- 

At the end of this section we would like to present an interesting observation. The 
formula analogous to Eq. (|25|) for the case of Qa has the form 



o {1) = NcM a (0) (26) 

a A, sea, lowest 3>2,1X^J "^> sea >l° wes t ' \ / 

Note that the coefficient between the 1/N C correction and the leading- term in Eq. 



is exactly twice the analogous coefficient in Eq. (|25|) . This is not a coincidence. In the chiral 
limit, — > 0, the sea contribution to a is divergent as l/m^ P JT^ , |26[| . This divergent 
piece is proportional to g\, and this result is a general, model independent feature of the 
chiral limit. Therefore, one has a sea = ^-g% where 7 is a numerical constant. Hence, if in a 
model one corrects g& by including 1/N C terms, analogous correction should come out in the 
calculation of 1/N C corrections to a. Consider the formal case of solitons with a very large 
size R. In that limit, qa consists only of the sea contribution, which in turn is saturated by 
the lowest term in the gradient expansion: qa = fiU,sea,iowest- Using Eq. (|26|) and the general 
considerations discussed above, we obtain in the double limit of — > and R —>■ 00 the 
following result 

_ 7 2 _ 7 / (0) , (l) \ 2 , nh /N 2\ 

C*sea,lowest — 9 A ~ _ [ S'A.sea.lowest + fl'A.sea.lowest j ' Ly l i /- /V cJ 

= "seagull + Seagull + 0(1/N%). (27) 

The first term in the last line of is our «sea lowest; the second term is ag^ lowest , with exactly 
the same factor as in Eq. (|25|) . This shows that our method of including rotational 1/N C 
corrections is consistent with chiral physics. 

Analogous consistency conditions have to hold for other observables which behave as 
fi^/m?,-, e.g. the isovector magnetic mean squared radius of the nucleon. 

Another digression concerns the structure of the rotational 1/N C corrections. Collective 
operators which arise in evaluating corrections to or the isovector magnetic moment have 
the same structure as collective operators entering at leading-A^ level. Formally, this is 
because [J a , Dl c ] = ie abd D\ c . Thus, ratios of these observables (also for non-diagonal matrix 
elements) for various members of the flavor SU(2) multiplet (i.e. for and A) are preserved 
at the 1/N C level, and obtain corrections starting at the level. This complies to the 



consistency condition of Ref. |39|. 
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V. EFFECTS OF THE A-N MASS SPLITTING 



Before comparing our results to data, we note that there is another very important effect 
which has not been included, namely the iV-A mass splitting. As discussed in Ref. |T2],|2T) 



in the limit N c — > oo our perturbative approach has the nucleon and the A degenerate in 
mass. At the hadronic level, the seagull term results from diagrams with intermediate pion- 
nucleon or pion-A states. Because of the N-A degeneracy, the contribution of the pion-A 
states is strongly overestimated in the hedgehog treatment. One may crudely include the 
effect of N-A mass splitting by reducing the seagull contribution by about 30% |I2|JSH]. A 
consistent treatment of this effect would require the use of a better collective quantization 
method, but this does not seem feasible for the NJL model. As shown in Section [III] more 
than 80% of the sea contribution cc is carried by the seagull term. Therefore, as suggested 



in Refs. [|T2"| , p6"f , we simply multiply the sea contribution by a factor of 2/3 to mimic the 
effect of N-A mass splitting. 

Since a£][ lowest is proportional to a seagu n, we expect that it obtains analogous corrections 
from A-N splitting as the leading- N c part. Therefore we also multiply it by a factor of 2/3 
before comparing to data. 

VI. FINAL RESULT AND CONCLUSION 

Our numerical results are collected in Figs. 1,3. They are plotted as a function of the 
constituent quark mass M. Figure 1 shows the sea and valence contributions to the leading- 
Ac nucleon electric polarizability a^ ^ as well as to the rotational 1/N C correction a^. The 
total polarizability is presented on Fig. 3 compared to the experimental value. The dashed 
line corresponds to a' = a^l + afj + lowest + a£j[. The full line shows the total a, crudely 
corrected for the effects of the N-A mass splitting by multiplying the sea contributions by a 
factor of 2/3 (Sec. |VJ) : a = §a£) + afj + §a2a,iowe S t + «S- F °r values of M in the physically 
relevant range 400-450 MeV the changes of a with M are small. Choosing for M = 420 MeV 
(value which gives good agreement for other nucleon observables |HJ|21J) we obtain for the 



total nucleon electric polarizability the value a = 19 x 10~ 4 fm 3 , with §0^, ck^, lowest) 
and ttyai contributing 36%, 29%, 17%, and 17% respectively. Our result is still too large 
compared to the experimental number a ex p = 9.6 ± 1.8 ± 2.2 x 10~ 4 fm 3 |0-|5|. Keeping 
in mind the approximate nature of the model and techniques used, our model prediction is 
within expectations. 

Note, that in our calculation, after including the rotational 1/N C effects, we have the 
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correct value of qa- This is not the case of other models [|,|T(],[T3|] , and apparent agreement 
of a with experiment in these models results from having too low g^. 
The main points of this paper can be summarized as follows: 

1) Dirac sea effects dominate the electric polarizability of the nucleon. In the gradient 
expansion the sea effects are dominated by the seagull contribution. This explains how this 
term arises in other models which do not have the explicit Dirac sea. 

2) We have calculated rotational 1/N C corrections, and found that they are important 
and of expected size. 

3) We have crudely estimated the N-A splitting effects. 

4) With all above effects included, numerical predictions of the NJL model are better 
than in other soliton models (note that we have the correct value of g/S) but are still about 
a factor of 2 too large compared to experiment. This indicates that the NJL soliton is too 
"soft" and other effects are still needed in order to describe the polarizability properly. 

5) Calculation of the magnetic polarizability, or the proton-neutron splitting of the elec- 
tric polarizability would require a much greater effort in the NJL model than in the cx-model 
calculation of Ref . |JT^JT3] . This is because of the necessity to subtract zero-modes. 
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FIGURE CAPTIONS 



Fig. 1. 

The electric polarizability of the nucleoli, a, as a function of the constituent quark mass M. 
Leading N c Dirac sea contribution a^l (dashed line); leading N c valence contribution a^j 
(dash-dotted line); rotational 1/N C correction to the sea contribution ctgga lowcst (dotted line) 
and to the valence one ct^j (dash- dot- dot). 

Fig. 2. 

The ratio of the gradient expanded electric nucleon polarizability to the full leading-iV c 
= at^l + cn^j for different soliton sizes R. The dash-dotted line represents the leading 
order in the gradient expansion a sea guii an d the dashed line the seagull plus the next term 
in the gradient expanded a. 

Fig. 3. 

The total electric polarizability of the nucleon, a, as a function of the constituent quark 
mass M. The dashed line shows the full a' = a^l + a^] + a£^ lowest + a^j and the full line 
a = |a^el + a lai + | a sea,ieadin g + a iab where the factor I in front of crudely corrects for 
the overestimation of the A contribution (Sec. |V|). 
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